SECURITY  CL.  ASSlFiCATiQN  0*  This  PAGE  <’P7t#o  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

i  report  number 

\Z.  GOVT  ACCESSION  MO.| 

3.  RECIPIENT'S  CATALOG  num»ER 

NRL  Memorandum  Report  4659 

i 

i  ... 

4  TITLE  fend  Subtitle) 

1 

S  TYPE  OF  REPORT  4  PERIOO  COVERED 

Interim  report  on  a  continuing 

PRELIMINARY  NUMERICAL  SIMULATION  OF  IR  STRUC-  ! 

NRL  problem. 

TURE  DEVELOPMENT  IN  A  HYPOTHETICAL  URANIUM 

ft.  PERFORMING  org.  report  number 

RELEASE 

7  AUTHOR^ 

ft.  CONTRACT  or  grant  number^ 

E.  Hyman*,  M.  Mulbrandon,  S.  L.  Ossakow,  and 

B.  E.  McDonald  *♦ 

t.  RERFORMING  ORGANIZATION  NAM  £  AND  AODRESS 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  S  WORK  UNIT  NUMBERS 

Naval  Research  Laboratory 

Washington,  DC  20375 

62074H;  47-0917-0-1 

ii-  controlling  office  name  ano  aooress 

! 

12.  REPORT  DATE 

Defense  Nuclear  Agency 

November  16, 1981 

Washington,  DC  20305 

13  NUMBER  OF  PAGES 

51 

i 


MONITORING  agency  NAME  4  AOORESSri/  dillarant  tram  Controlling  Otfice) 


>3.  SECURITY  CLASS,  tot  title  raport) 

UNCLASSIFIED 


DISTRIBUTION  STATEMENT  (o I  thim  Raport) 

Approved  for  public  release;  distribution  unlimited. 


^  i 

'  |  Y  -  £  Y 


IT.  DISTRIBUTION  STATEMENT  ( of  the  abstract  an  tat  ad  in  Block  20,  II  dillarant  hoot  Raport) 


IB.  SUPPLEMENTARY  NOTES 

♦Present  address:  Science  Applications,  Inc.,  McLean,  VA  22102 

♦♦Present  address:  Naval  Oceanographic  Research  and  Development  Activity,  Bay  St.  Louis,  MS  39520 
This  research  was  sponsored  by  the  Defense  Nuclear  Agency  under  Subtask  S99QAXHI,  work  unit 
00001,  and  work  unit  title,  “IR  Structure.” 


19-  KEY  WOROS  (Continue  on  ravaraa  mi  da  if  nacaaaery  and  Identify  by  block  numbar) 


IR  Structure 
Uranium  release 
Numerical  simulation 
Isophots 


Power  spectrum 
Parallax  effects 
PHARO  code 
LWIR 


20.  ABSTRACT  (Continue  on  reaerae  alda  II  nacaaaary  and  Identity  by  block  numbar) 
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altitude.  Of  interest  is  the  LWIR  emission  from  uranium  oxide  ions,  induced  by  sunlight  and  earthshine. 
Assuming  a  one-level  fluid  and  integrated  Pedersen  conductivity  ratio  (cloud  to  ambient)  of  5,  the  cloud 
diffuses  along  the  magnetic  field  in  time  and  rapidly  develops  gradient  drift  striations  transverse  to  the 
magnetic  Reid  and  across  the  direction  of  the  relative  velocity  of  neutrals  to  ions.  Using  the  PHARO 
code,  we  produce  contour  plots  of  radiance  in  the  11-14  u  region  to  simulate  observations  from  various 

(Continues) 


I  JAW  71  1473  COITION  Or  1  NOV  #S  IS  QSSOLCTC 

S/N  0102-014-  ««0  1 


SECURITY  CLASSIFICATION  OF  THIS  RAOE  felon  O.t.  InlwMJ 


SCCUWITV  CLASSIFICATION  or  TwK  r»*—v  Omtm  £*»•  r*4) 


20  ABSTRACT  (Continued) 

locations  both  along  and  transverse  to  the  magnetic  field.  Contours  of  radiance  vary  in  Intensity  and 
character  from  different  vantage  points.  Power  spectra  are  calculated  to  quantify  the  development  of 
spatial  irregularities  in  the  radiance.  We  find  parallax  to  be  an  important  effect  which  alters  both  the 
appearance  of  the  radiance  contours  and  the  spectral  characteristics.  This  inherently  3D  effect  could 
not  be  observed  in  previous  2D  simulations.  Altering  the  distance  from  observer  to  cloud  or 
changing  slightly  the  direction  of  observation  can  produce  drastic  changes  in  the  observations. 
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PRELIMINARY  NUMERICAL  SIMULATION  OF  ER  STRUCTURE 
DEVELOPMENT  IN  A  HYPOTHETICAL  URANIUM  RELEASE 


I.  INTRODUCTION 

The  significance  to  U.S.  defense  systems  of  long  wave 
infrared  (LWIR)  emissions  from  metallic  oxides  in  the  debris  from 
a  high  altitude  nuclear  explosion  (HAN'E)  is  an  important  unanswered 
question.  Of  most  concern  are  structured  or  striating  media  since 
they  degrade  optical  systems  at  intensities  many  orders  of  magnitude 
lower  than  smooth  media.  Those  debris  species  which  result  in  long 
lived  ions  are  of  most  interest  because  they  are  the  ones  most 
likely  to  structure.  For  this  reason  uranium  deserves  careful 
scrutiny.  Because  of  the  low  ionization  potential  of  UO ,  uranium 
reacts  with  atomic  oxygen  via  an  associative  ionization  process  to 
form  U0+.  Thus,  this  molecular  ion  is  stable  to  dissociative  re¬ 
combination,  and  deionization  requires  a  three  body  reaction. 

Uranium  deposited  in  the  high  altitude  atmosphere  (>100  km)  in  a 
nuclear  burst  will,  therefore,  induce  very  long  lived  ionization. 

The  emission  characteristics  of  UO+  have  not  been  measured, 
although  on  general  theoretical  grounds  a  very  large  line  density 
is  expected  in  the  ll-14p  region.  A  laboratory  experiment  to 
measure  the  U0+  spectrum  (wavelengths  and  band  strengths)  in  this 
region  of  the  Llv'IR  has  been  described  by  Bien.1-^  It  has  not  yet 
been  carried  out. 

Because  of  the  potential  importance  of  uranium  debris 
from  HAXE  events  in  the  performance  of  defense  optical  systems, 
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the  RAAE  Division  of  the  Defense  Nuclear  Agency  has  suggested 

that  numerical  simulations  be  performed  to  model  high  altitude 

uranium  releases.  Based  on  reasonable  estimates  of  oscillator 

strengths  for  transitions  in  the  LWIR,  these  simulations  will 

predict  structuring  characteristics  of  the  hypothetical  plasma 

cloud  released  in  the  ionosphere.  In  this  paper  we  will  present 

results  of  a  three  dimensional  simulation  of  a  uranium  release. 

2 

Using  the  PHARO  code  we  will  show  plots  of  radiance  (watts/cm 
ster)  in  the  ll-14y  region  as  a  function  of  time,  (assuming  the 
emission  to  be  optically  thin) ,  viewing  both  along  and  trans¬ 
verse  to  the  magnetic  field.  These  plots  will  show  the  develop¬ 
ment  of  structure  in  time.  In  addition,  we  will  show  the  power 
spectra  for  the  structure,  giving  the  distribution  of  spatial 
wavelengths  in  the  emitting  band.  Parallax  effects  will  also  be 
described. 

II.  HYPOTHETICAL  EXPERIMENT 

We  will  simulate  a  hypothetical  release  based  on  the 

parameters  presented  at  the  DNA  Infrared  Program  Meeting  at  the 

.  (2) 

Naval  Research  Laboratory  in  January  1980  by  W.  Reidy; v  A 

release  at  200  km  altitude  in  daylight  with  an  atomic  uranium 

24 

content  of  2  kgm  (^5  x  10  atoms)  is  modeled.  Initially  the 
uranium  vapor  will  expand  radially  with  a  thermal  velocity  some¬ 
what  less  than  1  km/sec.  The  expanding  uranium  vapor  would  be 
slowed  when  it  has  interacted  with  an  equal  mass  of  the  atmo¬ 
sphere  (^1.2  km).  However,  before  this  occurs  it  will  interact 

with  an  equal  number  of  atomic  oxygen  atoms  (>.75  km)  .  Assuming 

9  —  3 

an  oxygen  density  ^3  x  10  cm  and  a  reaction  rate  coefficient 
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of  4  x  10  ^cm'Vsec,  ^  the  ionization  will  occur  with  an  e- 
folding  time  of  ^.8  sec.  This  will  restrict  further  growth  of 
the  cloud  transverse  to  the  magnetic  field,  and  the  cloud  will 
initially  develop  a  sausage-like  shape  with  a  radial  dimension 
transverse  to  the  magnetic  field  ^.75  km.  However,  2x2  convec¬ 
tion  will  elongate  the  cloud  in  the  direction  of  the  relative 
velocity  of  the  neutral  atmosphere  to  the  ion  cloud,  and  will 


narrow  it  in  the  other  direction  transverse  to  the  magnetic 


field. 


(4) 


Because  of  the  presence  in  the  normal  atmosphere  at  200  km 

of  a  substantial  amount  of  molecular  oxygen,  which  may  not  be 

present  in  the  nuclear  disturbed  atmosphere,  a  further  reaction 

will  occur;  UO+  will  react  with  O2  to  form  UO*.  Assuming  a 

-9  3  (55 

reaction  coefficient  of  2  x  10  cm  /sec  and  an  0 2  density  of 

8*3  4“ 

v4  x  10  cm  the  e-folding  time  for  conversion  of  uo  to  U02  is 

just  1.2  sec.  Thus,  after  the  first  several  seconds  the  predomi¬ 
nant  ion  will  be  UO^.  This  will  have  little  effect  on  the  dynamics 
and  structuring  of  the  ion  cloud  but  the  LWIR  spectrum  will  be 
different,  of  course.  Because  in  the  nucleor  case  the  02 
molecules  may  be  dissociated,  the  transition  from  UO+  to  UC>2 
may  take  much  longer.  On  the  other  hand,  even  in  a  high  altitude 
burst  much  of  the  debris  will  typically  be  deposited  at  lower 
altitudes  (^100  km)  where  there  will  be  large  amounts  of  02 . 

Metal  oxides  are  excited  from  the  ground  state  to  higher 
electronic  levels  by  visible  radiation  from  sunlight.  These 
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excited  levels  rapidly  decay,  in  general,  to  excited  vibrational 

levels  of  the  ground  electronic  state  which  radiate  in  the  LWIR. 

In  addition,  there  is  direct  excitation  of  vibrational  levels 

caused  by  earthshine  at  LWIR  wavelengths,  leading  to  LWIR  emission 

as  these  levels  decay.  Typical  emission  rates  for  metal  oxides 

_  -> 

in  the  LWIR  are  'vlO  "  photons/molecule/sec  due  to  earthshine,  and 
^10  ^  photons/molecule/ sec  due  to  sunlight.  There  has  been 

some  speculation  that  the  emission  rates  for  U0+  and  UO*  may  be 
substantially  higher  than  these  values^*  because  of  the  presence 
of  low-lying  electronic  levels,  but  that  has  not  been  established. 
In  the  case  of  aluminum  oxide  the  rate  is  v.2  photons, molecule/sec 
both  for  sunlight  and  for  earthshine.  LWIR  emission  is  expected  to 
be  in  the  ll-14u  region  but  the  results  in  this  paper  do  not  depend 
on  this  specific  wavelength  interval. 

III.  THE  SIMULATION 

The  model  assumes  a  one  ionic  species,  one-level  fluid  and 

(  8  ) 

follows  the  analysis  used  for  the  Avefria  barium  cloud  releases. 

The  conductivity  ratio  (U0  +  plasma  cloud  to  ba  :kgi  ound  ionosphere 

Pedersen  conductivity  ratio),  as  defined  in  Reference  8,  is  taken 

to  be  5  which  is  based  on  an  assumption  of  301  efficiency  of 

vapori ration  of  the  uranium.  Briefly,  the  conductivity  ratio  5 

■>  4 

was  obtained  as  follows.  The  5  x  10"  atoms  of  U  was  assumed  to 
be  30 o  vaporized  and  the  resulting  U0+  was  gaussianly  distributed 
over  a  sphere,  with  a  gaussian  radius  of  0.“5  km.  Using  a  UO 
ion-neutral  collision- frequency  (.v^)  of  0.2  s  ^  and  a  gyro- 
frequency  (O.)  of  18. S  s  *,  the  integrated  Pedersen  conductivity 
of  the  UO*  cloud  was  found  to  be  50  mhos.  For  the  daytime  high 


4 


latitude  ionosphere,  a  value  of  6  mhos  for  the  integrated 
Pedersen  conductivity  was  taken. 

The  altitudes  and  sices  of  the  ionized  clouds  are  similar 
in  the  two  releases  as  are  the  striation  formation  mechanisms  and 
lead  to  the  rapid  development  of  gradient  drift  striations. 

Whereas,  in  the  Avefria  simulation,  however,  the  F  region  ionospheric 
plasma  cloud  was  reduced  to  two  dimensions  (transverse  to  the  magnetic 
field)  by  performing  f ie Id- 1 ine - integra t ions  over  the  Pedersen  con¬ 
ductivity,  here  in  the  model  this  is  not  done  and  the  cloud  also 
diffuses  in  time  along  the  magnetic  field,  providing  a  three  dimen¬ 
sional  simulation.  This  enables  one  to  view  the  cloud,  using  the 
f  9 1 

PHARO  code,  J  from  an  arbitrary  direction  and  to  automatically 

incorporate  any  effects  of  parallax  due  to  the  three-dimensionality 

of  the  cloud.  A  PHARO  plot  is  obtained  by  summing  the  volume 

emission  rates  along  an  array  of  rays  originating  at  the  camera 

location.  A  2D  array  of  radiance  values  is  determined ,  from  which 

a  contour  plot  of  the  radiance  is  constructed. 

Thus,  the  cloud  is  spherical  for  the  first  few  seconds 

until  it  becomes  ionized.  From  then  on,  after  this  time  t^, 

diffusion  transverse  to  the  magnetic  field  can  be  neglected. 

Along  the  magnetic  field,  we  assume  an  initial  density  distri- 

•>  *> 

bution  n(t^)  ^  exp(-z'/L"),  where  z  is  the  coordinate  along  the 
magnetic  field  (B)  direction  and  L  %  ."5  km.  With  a  parallel  to 
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B  diffusion  constant  K  =  c“/v .  ,  where  c  'v  .  3  km/sec  and 

—  s  in  s 

+  - 1  . 1  ~> 
v>in(UO  )  'v  .2  sec  ,  we  have  Kji'v  4.15  x  10  kmVsec.  The  density 

diffuses  in  time  according  to 


0 


(1J 


At  time  t  then,  we  will  have 


r  1  i 

1/^  /  ,2 

1 

|_  1  +  b  ( t  - 1  i  ) 

expj  -  — 

V  L" 

1  +  b  ( t  - 1  i ) 

with  b  =  4  K  /  L“  \  3. 1  sec' 1 . 

The  above  analysis  is  strictly  correct  only  for  a  uniform 
atmosphere.  For  a  release  in  equatorial  regions  the  z  direction 
is  horizontal  and  there  is  no  problem.  At  high  latitudes,  where 
the  magnetic  field  is  nearly  vertical,  diffusion  will  be  signifi¬ 
cantly  faster  in  the  upward  direction  than  in  the  downward  direction. 
\ow,  the  diffusion  scale  size  grows  like  the  square  root  of  the  time 
and  at  200  km  is  x  13.4  km  at  100  sec.  Therefore,  at  times  of  the 
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order  of  100  sec  this  effect  is  not  major  (atmospheric  scale  size 

at  200  km  %  45  km).  We  have,  for  this  reason,  not  corrected  for 

this  effect  in  our  simulation. 

The  above  analysis  is  performed  while  simultaneously 

following  the  perpendicular  motion  and  structuring  of  the  U0+ 

(  8  ) 

plasma  cloud  according  to  the  one-level  set  of  equations  • 


M-  •  •  t  -  iw > 


V  =  -  |  7  t  x  2 
—  B  x 


7 

X 


(S  Va3>) 


E 

— o 


7  £ 

x 


(5j 


where  £  ,  V,  K  ,  $>,  B  and  E  are  the  integrated  plasma  cloud 
i UO+ !  Pedersen  conductivity  (y  =  f  a  (x,y,z)dz,  where  o  is  the 
Pedersen  conductivity)  and  is  a  function  of  x,  y  (coordinates 
perpendicular  to  EH,  perpendicular  plasma  velocity,  relative 
to  the  ambient  perpendicular  plasma  drift,  the  cross-field 
(x  B)  diffusion  coefficient,  the  cloud  -  induced  potential,  magnetic 
field  strength,  and  ambient  perpendicular  electric  field  in  the 
rest  frame  of  the  uniform  neutral  flow,  respectively  and  the 
electrostatic  approximation  (E  =  )  has  been  made.  In 

addition  -  =  l  2kT/nu  )  iven+ve  )  where  k,  T,  nw  ,  vgn, 

.  e  and  are  Boltzmann's  constant,  the  plasma  temperature, 

ion  mass,  electron-neutral  collision  frequency,  electron-ion 
collision  frequency,  electron  gyro  frequency  and  ion  gyrof requency , 
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1 

1 


respectively.  Values  of  E  =  SmV/'m,  B  =  0.5  gauss  (and  resulting 

7 

in  a  cEq/B  velocity  of  100  m/s),  and  ^  i  1  m"/sec  were  chosen  for 
the  simulations.  The  computational  mesh  consisted  of  162  x  82 
grid  points  in  the  x  and  y  directions  (see  Figure  1) .  The  plasma 
density  distribution  n(.x,y)  at  each  time  was  unfolded  from  (3)  - 
[5  i  by  distributing  the  density  along  z  according  to  the  gaussian 
(with  a  time  dependent  gaussian  length)  given  by  (1)  and  (2),  using 
;^n(UO+i  and  £).(UO+),  integrating  over  z  and  setting  the  resultant 
integral  equal  to  £  (x , y ) . 

In  Figure  1  we  indicate  the  coordinate  system  fixed  in  the 
ionized  cloud,  and  illustrate  the  formation  of  striations  trans¬ 
verse  to  the  magnetic  field,  B,  and  across  the  direction  of  the 
relative  velocity  of  the  neutral  atmosphere,  v  ,  to  the  ionized 
plasma.  (In  this  system  E^  x  B  would  be  in  the  -  x  direction). 

The  calculat ional  mesh  follows  the  striating  edge  and  has  an 
outflow  boundary  condition  on  the  negative  x  boundary.  Thus,  in 
time  more  and  more  ions  are  lost  from  the  calculation  out  of  the 
"tail".  Figure  1  also  exhibits  the  locations  of  3  potential 
observation  directions  for  infrared  detectors  (cameras,).  In  the 
next  section  we  will  show  plots  of  ir  emission  contours  as 
observed  from  each  of  the  camera  locations.  The  emission  rate 
is  assumed  to  be  0.11  photons/ ion/sec  in  the  11 -lip  band. 

IV.  RESULTS 

a.  PHARC  Plots 

*> 

Figures  2-8  show  the  contours  of  radiance  in  watts/cm“/ster 
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at  several  times  as  observed  from  camera  1,  placed  a  distance 
away  of  200  km.  Camera  1  looks  along  the  magnetic  field. 

Figure  2  shows  a  PHARO  plot  of  the  radiance  50  sec  after  the 
release.  The  abscissa  and  ordinate  are  in  degrees  and  are  re¬ 
spectively  along  the  x  and  y  directions  in  the  coordinate  system 
of  Figure  1.  An  82  x  82  mesh  has  been  used  in  the  PHARO  part  of 
the  calculation  in  the  camera  plane  of  view.  The  contour  values, 
defined  in  the  headings  of  the  plot,  potentially  range  from 

1.125  x  10  watts/cm-  ster,  (the  "9”  contour)  to  0.0  (the  "0" 

-  8  7 

contour)  in  steps  of  1.25  x  10  watts/cm"  ster.  The  largest 

contour  appearing  in  Figure  2  is  the  "8"  contour  which  represents 

_  •*> 

a  radiance  of  1.00  x  10  '  watts/cm-  ster.  A  key  at  the  bottom  of 

■> 

the  plot  gives  the  value  of  some  of  the  contours  in  watts/cm- 
ster.  At  60  sec  there  is  no  sign  of  the  fingering  behavior  which 
becomes  so  evident  at  later  times.  In  Figure  3,  80  sec,  the 
striating  behavior  has  clearly  begun  developing.  The  structuring 
intensifies  at  100,  120,  140,  160  and  180  sec  shown  in  Figures  4, 

5,  6,  _  and  8  respectively.  The  overall  radiance  has  not  diminished 
with  time. 

Camera  2  looks  along  the  relative  velocity  vector  of  the 
neutral  atmosphere  and  the  ion  cloud  from  a  distance  of  200  km. 

In  Figure  9  we  exhibit  the  PHARO  plot  at  140  sec  as  observed  from 
camera  2  with  a  sufficiently  large  field  of  view  that  the  entire 
cloud  is  in  view.  While  this  plot  is  useful  to  exhibit  the  entire 
cloud,  there  is  not  enough  resolution  to  distinguish  the  contours. 
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For  this,  we  return  to  the  high  resolution  in  Figure  10,  which 

is  at  bO  sec.  Only  a  portion  of  the  cloud  is  in  view  here,  of 

course.  The  magnetic  field  in  Figure  10  is  along  the  ordinate, 

and  the  abscissa  is  the  v  coordinate  of  Figure  1.  Here  the  maximum 

"8  0 

contour  value  observed  is  1.0  x  10  watts/cm"  ster.  The  values  of 
radiance  determined  from  camera  2  are  substantially  lower  than  those 
from  camera  1  since  from  this  vantage  point  lines  of  sight  do  not 
pass  through  as  much  of  the  cloud  as  in  the  end  view  observed  by 
camera  1.  Furthermore,  from  camera  2  we  are  looking  up  the  "tail", 
along  the  x-direction,  and,  as  we  pointed  out,  some  of  the  ions 
have  left  the  calculat ional  mesh.  This  means  that  not  only  are 
the  radiance  levels  too  small  as  observed  from  this  direction,  but 
any  nonuniformity  in  the  "tail"  has  been  lost.  On  the  other  hand, 
there  is  no  fine  structure  lost,  as  this  is  all  in  the  leading  edge. 
Our  calculations  indicate  that  approximately  1/3  of  the  ions  have 
left  the  mesh  by  140  sec. 

Note  the  near  perfect  symmetry  at  bO  sec,  corresponding  to 
the  symmetry  observed  from  camera  1  at  this  time.  In  Figures  11-14 
we  observe  the  changes  in  time  from  the  camera  2  vantage  point.  At 
100  sec.  Figure  11,  we  see  the  maximum  contour  begin  to  shift  to 
the  left,  corresponding  to  the  maximum  in  Figure  4  shifting  down¬ 
ward.  A  relative  minimum  forms  just  to  the  right  of  the  central 
vertical  axis  corresponding  to  the  depression  in  Figure  4  just 
above  the  central  horizontal  axis.  These  effects  become  more 
pronounced  at  succeeding  times.  Particularly  interesting  from 


this  direction  is  the  alternation  in  intensities  as  we  scan 
horizontally  across  the  field  of  view  (Figures  12-14).  There 
is  a  gradual  reduction  with  time  in  the  radiance,  also,  as  more 
of  the  cloud  diffuses  along  the  magnetic  field  out  of  the  field 
of  view,  and,  also,  out  of  the  calculat  iona 1  mesh  (negative  x 
boundary).  The  former  is  a  real  physical  effect,  the  latter 
an  artifact  of  the  calculation. 

Camera  3  looks  in  the  direction  transverse  to  the  magnetic 
field  and  to  the  drift  direction  from  a  distance  of  200  km.  Again, 
in  these  views  (Figures  15-19)  not  all  of  the  cloud  is  included. 

The  ordinate  is  along  the  magnetic  field  (;)  and  the  abscissa  is 
along  the  relative  velocity  (x) .  The  steep  gradients  at  the  left 
of  Figure  15  correspond  to  those  on  the  right  in  Figure  2.  As 
from  camera  2,  the  radiance  decreases  with  time  as  the  cloud 
diffuses  out  of  the  field  of  view  but  there  is  no  loss  from  missing 
ions,  as  the  "tail"  is  not  in  the  field  of  view  of  camera  5.  Here 
the  structuring  leads  to  an  appearance  significantly  distinct  from 
either  of  the  other  two  views. 

b.  Power  Spectra 

Spatial  irregularity  power  spectra  have  been  produced  for 
each  camera  view.  For  example,  consider  the  radiance  as  viewed 
from  camera  1  at  a  distance  of  200  km  at  140  sec  (Figure  6). 

Using  a  mesh  dimension  of  32  x  82  we  construct  a  two  dimensional 
(2D)  power  spectrum.  This  is  converted  to  two  one  dimensional 
spectra:  the  x  power  spectrum  is  a  result  of  integrating  over  the 
v  components  of  the  wavenumber  (k.)  of  the  2D  power  spectrum;  the 
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>'  power  spectrum  comes  from  integrating  over  the  x  components 
(k  1.  Typical  results  are  shown  in  Figures  20  and  21.  The  power 
scale  on  the  ordinate  is  in  units  of  the  logarithm  of  the  square 
of  the  radiance.  The  abscissa  has  two  scales:  (lj  the  logarithm 
of  the  spatial  wavelength  in  km,  and  (.2)  the  logarithm  of  the  mode 
number,  k/k^,  where  k^  =  2tt/D.  In  the  y  direction  D  is  the  system 
size  =  2. "km.  By  the  system  size  we  mean  the  distance  across  the 
field  of  view  based  on  the  distance  between  the  camera  and  the 
cloud,  and  neglecting  the  dimensions  of  the  cloud.  In  the  x 
direction  D  is  twice  the  system  size,  because  of  the  Neumann 
boundary  condition  in  x. 

A  typical  power  spectrum  plot  can  be  conveniently  divided 
into  two  segments.  The  low  mode  number  region  at  the  left  contains 
most  of  the  power  and  characterizes  the  gross  structure.  The 
positions  of  peaks  in  this  region  are  a  measure  of  the  size  of  the 
structure,  e.g.  in  Figure  21  there  is  a  clear  peak  at  a  wavelength 
just  under  0.5  km.  This  characterizes  the  typical  widths  of  the 
fingers  in  Figure  0.  The  second  region  is  the  asymptotic  regime, 
which  begins,  roughly  speaking  at  a  wavelength  corresponding  to 
the  smallest  scale  size  of  the  gross  structure.  This  region  can 
be  fit  approximately  by  a  straight  line  (in  the  log  of  the  power 
vs.  the  log  of  the  mode  number).  The  slope  of  that  fit,  X ,  is 
the  asymptotic  slope  of  the  power  spectrum.  The  value  of  the 
slope  is  determined  by  the  steepest  gradient  in  the  field  of  view. 

A  value  of  -2  corresponds  to  a  square  wave;  steeper  slopes  to  more 
gradual  gradients.  We  have  discarded  the  highest  15*  of  the  modes 
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from  this  calculation  since  the  highest  modes  are  subject  to 
severe  'aliasing'  effects  and  to  computational  errors. 
resolution  we  used  (82  x  82)  was  dictated  by  the  cost  of  the 
PHARO  calculations  and  should  be  considered  only  marginally 
adequate.  Thus,  the  values  of  ,\  should  be  considered  illustrative 
but  not  definitive.  Figures  20  and  21  do  indicate  steeper  gradients 
in  the  x  direction  than  in  the  y  direction. 

We  find  that  the  steep  gradients  are  established  fairly 
early,  certainly  by  100  sec.  However,  as  time  proceeds  the 
number  of  steep  gradients  increases.  This  shows  up  as  an  increase 
in  the  power  level  of  this  part  of  the  spectrum,  but  with  no  change 
in  slope.  Since,  as  we  have  said,  very  little  of  the  power  is  in 
the  asymptotic  part  of  the  spectrum  there  need  be  very  little 
change  in  the  power  spectrum  at  the  lower  modes  to  accommodate 
this  increase. 

From  camera  2  we  view  along  the  x  axis.  Thus,  the  only 
power  spectrum  that  is  meaningful  is  the  y  power  spectrum,  in 
which  the  z  components  have  been  integrated  over.  The  ion  dis¬ 
tribution  in  the  z  direction  is  strictly  gaussian,  and,  thus, 
uninteresting.  Figure  22  shows  the  y  power  spectrum  from  camera 
2  at  140  sec  at  a  distance  of  200  km.  Refer  to  Figure  13.  While 
there  is  clearly  mode  structure  here  in  the  asymptotic  region  the 
value  of  N  is  not  too  different  from  that  in  Figure  21. 

From  camera  3  we  view  along  the  y  axis.  Only  the  x  power 
spectrum,  in  which  the  z  components  have  been  integrated  over,  is 
of  interest.  Figure  23  shows  the  x  power  spectrum  from  camera  5 
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at  140  sec  at  a  distance  of  200  km.  Refer  to  Figure  18.  Again, 
there  is  a  great  deal  of  structure  in  the  asymptotic  region  but 
the  general  falloff,  characte r i ted  by  th<  N  value,  is  not  too 
different  from  that  in  Figure  20. 

c.  Parallax  Effects 

Previous  simulations  ‘  have  been  two  dimensional  so  that 
parallax  effects  were  unobservable.  We  find  that  including  the 
thi rd • dimens  ion  has  a  profound  effect  both  on  the  radiance  con¬ 
tours  and  on  the  irregularity  spectral  characteristics.  We  will 
illustrate  this  in  two  ways,  by  showing:  (1)  how  these  properties 
vary  with  distance  and  (2)  how  they  vary  with  change  in  pointing 
direction. 

Consider  again  camera  1  at  140  sec  at  a  distance  of  200  km, 
Figure  6,  and  compare  with  Figures  24  and  25,  which  are  views  from 
camera  1  at  140  sec  at  respectively  50  km  and  600  km  from  the 
cloud.  Clearly,  at  50  km  parallax  has  spread  the  contours  greatly 
and  substantially  reduced  the  gradients.  At  60i  km  the  reverse 
has  occurred  and  the  gradients  are  larger.  This  is  borne  out  by 
the  power  spectral  results.  Figures  2b  and  27  are  the  x  power 
spectra  from,  respectively,  50  km  and  600  km  and  should  be  com¬ 
pared  with  Figure  20  for  200  km.  The  spectral  index,  N,  at  50  km 
is  -5.44,  at  200  km  is  -2.90,  and  at  600  km  is  -2.39  documenting 
the  clear  sharpening  of  the  steepest  gradients  as  we  recede  from 
the  cloud.  Figures  28  and  29  are  the  v  power  spectra  from, 
respectively,  50  km  and  600  km  which  should  be  compared  with 
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Figure  21  for  200  km.  The  v  spectral  index,  N ,  is  the  same  at 
50  km  and  200  km,  within  the  accuracy  of  the  fitting  technique, 
but  has  decreased  by  ^  10"  at  600  km.  This  is  because  in  the 
y  direction  there  are  steep  gradients  near  the  center  of  the  field 
of  view,  where  parallax  effects  are  small,  so  that  the  steep 
gradients  have  been  somewhat  preserved.  However,  if  we  compare 
the  power  levels  in  the  asymptotic  regime  at  50  km  and  600  km, 
at  some  intermediate  fixed  wavenumber,  we  find  the  power  down  by 
about  half  a  decade  both  in  the  x  and  y  spectra  (compare  Figures 
26  and  27  and  Figures  28  and  29)  at  50  km  relative  to  600  km. 

This  is  the  best  measure  of  the  effect  of  parallax  on  the  meas¬ 
ured  gradients. 

Suppose  now  we  imagine  that  camera  1,  at  140  sec  and  a 
distance  of  200  km,  is  rotated  about  the  x  axis  by  just  1°, 
toward  camera  3.  How  does  this  affect  the  radiance  contours? 
Figure  30  shows  the  resulting  contours,  which  have  been  greatly 
spread  out.  Although  the  maximum  radiance  values  have  been 
reduced,  compared  to  Figure  6,  the  irradiance  at  the  camera  is 
unchanged.  As  is  obvious  in  Figure  30,  the  greatest  change  has 
occurred  with  the  gradients  in  the  y  direction.  The  spectral 
index  of  the  y  power  spectrum  has  increased  from  %-4  to  'v-~. 

Figure  51  shows  the  effect  of  rotating  camera  1,  about 
the  y  axis  by  1°,  toward  the  x  axis.  Here  the  changes  from 
Figure  b  are  moderate.  However,  the  x  power  spectrum  index  does 
change  from  v-2.9  to  v-5.5  and  the  v  index  increases  by  about  20 
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Thus ,  when  we  rotate  toward  the  y  direction  we  lose  the  y  gradients 
and  when  we  rotate  toward  the  x  direction  we  lose  the  x  gradients. 

V.  SUMMARY 

The  severe  degradation  of  radio  signals  due  to  structuring 
of  electron  densities  in  the  atmosphere,  following  a  HANE  event, 
has  been  known  for  some  years  now.  A  great  deal  of  analysis  has 
been  performed  at  the  Naval  Research  Laboratory  (NRL)  to  under¬ 
stand  the  plasma  instability  mechanisms  that  underlie  this 
phenomenon,  as  well  as  related  phenomena  in  the  natural  ionosphere, 
such  as  spread  F.  It  is  clear  that  these  mechanisms  also  will 
result  in  structuring  in  the  emission  of  LWIR  radiation  and  con¬ 
sequently  is  potentially  of  serious  concern  to  the  U.S.  defense 
community.  Ultimately  we  want  to  use  our  theoretical  understanding 
of  striation  phenomena  along  with  our  capabilities  to  simulate 
HANE  phenomena  using  3D  magnetohydrodynamic  computer  models,  to 
permit  predictions  of  structuring  following  these  events. 

At  present  there  is  at  issue  the  question  of  what,  if  any, 
sources  of  LWIR  emission  are  likely  to  be  most  important  in  HANE 
events.  The  most  likely  candidates  at  present  are  Cl)  atomic 
oxygen  and/or  nitrogen  recombination  radiation  and  (2)  LWIR  emis¬ 
sion  from  metallic  oxides,  such  as  uranium.  To  shed  light  on  some 
of  the  unknowns  in  the  second  category  DNA  requires  a  reliable  model 
of  a  hypothetical  uranium  release  at  high  altitudes  to  clarify  the  emis¬ 
sion  wavelengths,  intensities,  and  structuring.  In  this  paper  we 
have  reported  on  the  results  of  a  numerical  simulation  of  a  uranium 
release  at  200  km  in  daylight,  making  reasonable  assumptions  about 
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the  values  for  some  of  the  unknown  parameters.  We  have  presented 
contour  plots  of  radiance  in  the  ll-14u  region  of  the  LWIR  as  it 
would  be  observed  from  detectors  placed  at  various  locations  along 
and  transverse  to  the  magnetic  field  and  at  various  distances  from 
the  cloud. 

This  enables  us  to  observe  the  development  in  time  of  the 
radiance  from  the  cloud  as  it  would  appear  from  any  vantage  point 
of  interest.  We  see  clearly  the  steepening  of  the  cloud,  the 
development  of  fingering,  and  how  this  is  related  to  the  power 
spectral  characteristics  that  are  observed.  Most  importantly, 
we  have  shown  that  the  three  dimensional  nature  of  the  cloud  has 
a  profound  effect  on  the  radiance  contours  and  on  the  observed 
spectral  characteristics.  Being  stationed  closer  or  observing 
slightly  off  axis  can  greatly  reduce  the  gradients  in  the  radiance 
of  the  cloud  and  change  the  power  spectral  characteristics  accord¬ 
ingly.  The  results  of  our  simulation  should  provide  a  theoretical 
understanding  of  a  hypothetical  experimental  release.  It  should 
also  be  recognized  that  the  simulations  we  have  presented  are  a 
first  cut,  i.e. ,  of  a  preliminary  nature.  Clearly,  improvements 
in  tn«=  numerical  simulations  presented  here  will  be  made.  Never¬ 
theless,  tnese  simulations  provide  some  interesting  results  on 
structuring  in  the  LWIR  regime. 
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Fig.  1  -  Schematic  view  of  the  uranium  cloud  and  the  coordinate  system, 
showing  diffusion  along  the  magnetic  field,  and  striation  formation  across 
the  direction  of  the  neutral-ion  relative  velocity  direction  v  .  Locations 
of  the  cameras  discussed  in  the  text  are  indicated. 
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Fig.  21  -  y  Power  spectrum  from  camera  1  (Figure  6  field  of  view) 
Refer  to  text  for  explanation  of  scales. 
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Fig.  22  -  y  Power  spectrum  from  camera  2  (Figure  13  field  of  view). 
Refer  to  text  for  explanation  of  scaxes. 


Fig.  23  -  x  Power  spectrum  from  camera  3  (Figure  18  field  of  view). 
Refer  to  text  for  explanation  of  scales. 
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Fig.  26  -  x  Power  spectrum  from  camera  1  (Figure  24  field  of  view). 
Refer  to  text  for  explanation  of  scales. 
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Fig.  27  -  x  Power  spectrum  from  camera  1  (Figure  25  field  of  view). 
Refer  to  text  for  explanation  of  scales. 
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Fig.  31  -  Camera  1,  rotated  1°  about  y  axis  toward  camera  2, 
radiance  contours  (watts/cm^ster)  200  km  from  the  cloud  origin, 
140  sec  after  release.  Abscissa  and  ordinate  are  in  degrees. 
Refer  to  text  for  contour  label  values  or  key  at  bottom. 
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OICY  ATTN  F.  FEHSENFELD 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 

P.  0.  BOX  983 
BERKELEY,  CA  94701 

OICY  ATTN  J.  WORKMAN 


KAMAN  TEMPO 

816  STATE  STREET  (P.O.  DRAWER  00) 
SANTA  BARBARA ,  CA  93102 
OICY  ATTN  DASIAC 
OICY  ATTN  TLM.  STEPHANS 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMAP^- 
01  CY  ATTN  MACK  STANTON 

GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.  0.  Box  6770 
SANTA  BARBARA,  CA  93111 
OICY  ATTN  JOHN  ISE,JR. 

OICY  ATTN  JOEL  GARBARINO 


INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY-NAVY  DRIVE 
ARLINGTON,  VA  22202 

OICY  ATTN  ERNEST  BAUER 
OICY  ATTN  HANS  WOLFHARD 

R  &  D  ASSOCIATES 

P.  0.  BOX  9695 

MARINA  DEL  REY,  CA  92091 

OICY  ATTN  FORREST  GILMORE 
OICY  ATTN  BRYAN  GABBARD 
OICY  ATTN  ROBERT  F.  LELEVIER 
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SCIENCE  APPLICATIONS,  INC. 

P.  0.  BOX  2351 
LA  JOLLA,  CA  92038 

01CY  ATTN  DANIEL  A.  HAMLIN 

MCDONNELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 

01CY  ATTN  R.  W.  HALPRIN 

MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  93101 

01 CY  ATTN  D.  SAPPENFIELD 
01CY  ATTN  R.  BOGUSCH 
01CY  ATTN  R.  HENDRICK 
01CY  ATTN  RALPH  KILB 
01CY  AITN  DAVE  SOWLE 
01CY  ATTN  F.  FAJEN 
01CY  ATTN  M.  SCHEIBE 

PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

01CY  ATTN  IRVING  L.  KOFSKY 

TECHNOLOGY  INTERNATIONAL  CORP. 

75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

C1CY  ATTN  W.  ?.  30QUIST 

SRI  INTERNATIONAL 
333  RAVENS WOOD  AVENUE 
MENLO  PARK,  CA  94025 

01CY  ATTN  WALTER  G.  CHESNUT 
01CY  AITN  CHARLES  L.  RINO 
01CY  ATTN  R.  HAKE,  JR. 

01CY  ATTN  D.  MCDANIEL 

UNITED  TECHNOLOGIES  RESEARCH  CENTER 
EAST  HARTFORD  CT  06108 

01CY  ATTN  H.  MICHAELS 

AERODYNE 

BEDFORD  RESEARCH  PARK 
EEDFORD,  MA  01730 

01CY  ATTN  C.  KOLB 

STEWART  RADIANCE  LABORATORY 
UTAH  STATE  UNIVERSITY 
1  DEANGELO  DRIVE 
BEDFORD,  MA  01730 

01CY  ATTN  J.  ULWICK 


VISIDYNE,  INC. 

5  CORPORATE  PLACE 
SO.  BEDFORD  STREET 
BURLINGTON,  MA  01803 

01CY  ATTN  W.  REIDY 
01CY  ATTN  J.  CARPENTER 
01CY  ATTN  C.  HUMPHREY 
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